Unambiguous discrimination between two unknown qudit states 
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We consider the unambiguous discrimination between two unknown qudit states in n-dimensional 
(n ^ 2) Hilbert space. By equivalence of unknown pure states to known mixed states and with 
the Jordan-basis method, we demonstrate that the optimal success probability of the discrimination 
between two unknown states is independent of the dimension n. We also give a scheme for a physi- 
cal implementation of the programmable state discriminator that can unambiguously discriminate 
between two unknown states with optimal probability of success. 
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I. INTRODUCTION 

As a recent development, the possibility of unambigu- 
ous discrimination between unknown quantum states can 
be potentially useful for many application in quantum 
communication and quantum computing. A universal de- 
vice that can unambiguously discriminate between two 
unknown states has been constructed by Bergou and 
Hillery [1]. It has three registers, labeled A, B and C, 
and each register can store a qubit that is in some arbi- 
trary state. In their work, it is assumed that register A 
is prepared in the state {ipi), register C is prepared in the 
states IV'2), and register B is guaranteed to be prepared 
in either jV'i) or \ip2)- Here, \ipi) and \ip2) are the states 
to be distinguished, which are both unknown 



iVi) =a|0) -1-611), |^.2> =c|0)+d|l), 



(1) 



where all the parameters a,6,c and d are arbitrary un- 
known complex variables satisfying the normalization 
equation |ap + |5p = 1 and |cp -I- |dp = 1. Furthermore, 
it is assumed that register B is prepared in the state jV'i) 
with probability rji and in the state 1-02) with probabil- 
ity r]2, such that 771 + 772 = 1, which guarantees that the 
state in register B is always one of these two states. The 
states viewed as a program which arc sent into registers 
A and C (called program registers) are called program 
states, while the unknown state to be conformed which 
is sent into the register B (called data register) is called 
data state. The device constructed here can measure the 
total input states 



1*1) = \^i)a\^i)b\^2)c, 
1*2) = Hi)a\'4'2)b\^2)c, 



(2) 



which are prepared with apriori probabilities 771 and 7/2. 
With the symmetry properties of the input states, this 
device will then, with some probabilities of success, tell us 
whether the unknown state in the data register B matches 
the state stored in program register A or C. 

In later works, a series of new devices have been in- 
troduced and widely discussed by other authors [2-7]. 



In these schemes we may have ha and nc copies of the 
states in the program register A and C respectively, and 
ub copies of the states in the data register B, then our 
task is to distinguish the two states 



1*1) = l^l)r"lV'l)r"l^2)g"^ 
i*2> = i^i)r"iV'2>r-i^2)r". 



(3) 



with the minimum error or with the minimum inconclu- 
sive probability, if the minimum error strategy or opti- 
mum unambiguous strategy is applied, respectively [8]. 

In both previous papers [2] and [8] , it has been proved 
that the optimal success probability of discrimination be- 
tween two unknown qubits is an increasing function of 
TLB, the number of copies in data register B. In this pa- 
per, we study the qudit states in n-dimensional Hilbert 
space and demonstrate that the optimal success proba- 
bility of discrimination between two unknown states is 
independent of the dimension n. Unlike the discrimina- 
tion between two known states, we cannot consider the 
subspacc spanned by the two states only, but we should 
consider the full n-dimensional space, as the two states 
are completely unknown to us. We adopt the equiva- 
lence between the discrimination of unknown pure states 
and that of known mixed states, and then reduce the 
problem to the unambiguous discrimination between two 
pure states with Jordan-basis method. Finally, we get 
the optimal success probability of the unambiguous dis- 
crimination between the two mixed states and give the 
detection operators. The detection operators are also 
applicable to the discrimination between the pure states 
without averaging over them. 

The organization of the paper is as follows. In Sec. II, 
we give a brief description of the programmable states 
discrimination, and adopt the equivalence between the 
discrimination of unknown states to that of known mixed 
states. In Sec. Ill, we introduced the Jordan basis for 
the mean input states. In Sec. IV, we find the optimal 
unambiguous discrimination and introduce its implemen- 
tation. Finally, a brief summary is given in Sec. V. 
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II. THE PROGRAMMABLE DISCRIMINATION 

In previous works, states such as \ip) = cos(0/2)|O) + 
sin(0/2)|l) have been considered [9]. These states he in 
a two-dimensional Hilbert space. Here we introduce a 
more gcnerahzed state 



(4) 



where is an arbitrary unknown complex variable. All 
these states span a n-dimensional Hilbert space "H, for 
which |i) form a mutually orthogonal basis. The two 
states we want to distinguish arc denoted by and 
\ip2)- As we mentioned before, a copy of each of the two 
unknown states is provided for the program registers A 
and C, denoted as \4>i)a and |'02)cj respectively. The 
state to be conformed is provided for the data register B 
as input. We assume that the state in the data register 
B is guaranteed to be prepared in \ipi) with probability 
771 and in 1^2) with probability = 1 ^ ?/i respectively. 
Thus we have two possible inputs 

l^-i) = |Vl)yl|V'l)B|V'2)c, 

1*2) = |V'1>^IV'2)b|V'2)c. (5) 

Here, and \ip2) completely arbitrary states 

n 

i=l 
n 

m = (6) 

i=l 

that remain unknown throughout the entire discrimina- 
tion process. 

Since the inputs and |^'2) arc unknown, and the 
states \ipi) and \ip2) making up these two inputs can 
change from preparation to preparation, it is only the 
pattern, middle states (B) matching the first (A) or the 
last (C) state that is preserved from one preparation to 
the next. Therefore, we can introduce the corresponding 
density operators 



Pl = {\^l){^l\}o 
P2 = {|*2)(*2|}o 



(7) 



where the average is taken over the entire parameter 
space of states 1-01) and |'02). Here, wc should notice that 
the optimal strategy for discrimination between the two 
states is the strategy that is optimal on average. That's 
to say, we can unambiguously discriminate between \ipi) 
and 1 0^2) as soon as we can unambiguously discriminate 
between pi and p2- 

To this end, we now define the space and operators 
that we will need. Let T-L be the n-dimcnsional space for 
the unknown states and then the full space is ®'^T-L. The 
space of symmetric states in H(E)'H is denoted by S, which 



is an 7T,(r7,-f l)/2-dimensional subspace. I'i'i) is an element 
of 5*1 = E (K) "H and |*2> is an element of ^2 = "H (8) S. 
Therefore, So = 5i n ^2, the intersection of Si and 5*2, 
is the space of symmetric states in the full space ^^H. 

50 is a subspace of dimension n{n + l){n + 2)/6, while 

51 and S'2 both n^(n -f- 1)/2. Let 5*3 be the subspace of 
(E)^H generated by 5*1 and S'2, and the dimension of 5*3 is 
ri(n-|-l)(5n— 2)/6. Let 5*4 be the orthogonal complement 
of So in Si , let S5 be the orthogonal complement of 

in ^2, and let 5*6 be the orthogonal complement of So in 
53. 

Clearly, the average in pi uniformly fills the symmet- 
ric subspace of A and B and the entire subspace of C, 
whereas the average in p2 uniformly fills the symmetric 
subspace of B and C and the entire subspace of A. There- 
fore, the corresponding density operators averaged over 
the unknown states, can be expressed as 

2 1 

Pi = — r-TT-Ps ® -Ic 

n[n + 1) n 

2 

11^ [n + i) 
1 



P2 =■ -Ia 



-Py. 



n n[n + 1) 

2f ^^iJ a ® Py, 



(8) 



where Ps is the projection onto S, and / onto T-L. Ps 
can be expressed as 



(9) 



(2\ 

where \u\j ) [i ^ j ~ l,...,n) arc the unique unit vectors 
in the symmetric subspace E, 



,(2)\ 



,(2)\ 



111) 
1 

71 



(|12) + |21)), 



(10) 



Consequently, we reduced the problem to discrimi- 
nating between the [n^(n -I- l)/2]-dimensional spaces Si 
and S'2 in S'3, which is equivalent to discriminating be- 
tween subspaces S4 and S5 in the [2n{n + l)(n — l)/3]- 
dimensional space 5*6. 



III. JORDAN BASIS 

It has been shown in [10] that two nonorthogonal sub- 
spaces T-Li and H2 of a Hilbert space can be unambigu- 
ously discriminated if we can find their canonical or Jor- 
dan bases [11]. The definition of the Jordan bases is as 



3 



follows. The set of orthogonal and normalized basis vec- 
tors {Iffi), in Hi and |/i„)} in 7^2 form 
Jordan bases if and only if 

{gi\hj) = % cos6'i, (11) 

where 9i arc the Jordan angles {Oi ^ 62 ^ ■■■ ^ On)- In 
the Jordan form, the full Hilbert space is decomposed 
into n orthogonal subspaces. Here, we denote the zth 
subspacc spanned by \gi) and \hi) by T^. Clearly, if the 
dimension of Ti is 1 , we can get cos 6i = 1 , and no 
further discrimination is possible in Tj. On the other 
hand, if cos^i < 1, \gi) and \hi) are linearly independent, 

I 



Then, the structure of the two density operators in Eq. 
(8); in particular the decomposition on the right-hand 

side, suggests that we consider and |a)(8)|u-j- ), 

where a = 1, n and 1 ^ i ^ j ^ n. These vectors form 
orthogonal bases for and S2, respectively. Thus the 



and they can be distinguished as two known pure states, 
which is familiar to us. Thus our task is reduced to some 
separate discriminations in each subspace, between two 
known pure states. 

Now, we will choose some bases to construct Jor- 
dan bases for the density operators pi and p2- First, 
we should notice that the fully symmetric subspace of 
the three states So must be common to both inputs. 
Here, we denote the basis for the subspace So by 
where l^i^j^fc^nis satisfied, and there are 
n{n -|- l)(n -|- 2)/6 these unique unit vectors, which can 
be expressed as follows, 



(12) 



unique unit vector in So can be expressed in terms 

of either Si or S2, since it is in both spaces. Here, we 
neglect the vectors |it,-jj,) that satisfy i ^ j ~ k, and a 
direct calculation shows that 



r m) 



1^15) = <! 



\ijk) + \ikj\) + \kij) 
^/3 

\ijk) + \]ik) + \]ki) 
V3 



if i = j = fc 

if i = J 7^ fc 



il i ^ j = k 



\ijk) + \]ik) + \ikj) + \kij) + \jki) + \kji) . 



V6 



\i i ^ j ^ k 



and 



iwf^)\3) + \llm ifz^jvfc 

= < ^lWfj')\k) + ^\jki) ifi^j^k , 



1K)I«S> 



l|fczj) 



if z = j 7^ fc 
ii i ^ j ~ k 

I 



l|zjfc) 

,(2)\ 



We can now introduce the vectors 



,(2) 



k)\J) 



\zjk) 
\kij) 



for i = j ^ k, 

(15) \g,,k) 

(16) \h,,k) 



jki) 
ijk) 



(13) 



(14) 



(17) 
(18) 
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for i ^ j = k, 



and 



\9ijk) 
\9ijk) 

\hi-jk) 
\Kjk) 



3-^3, (2),,,, 3 + 73 

"i7 a 



6 ' 
3-V3 (2), I . 

— I — I^V>b 



3 + %/3 



(19) 



3- \/3 



|fc>l' 



,(2)\ 



3 + V3, 
6 



,(2)\ 



6 

3- \/3 (2), 3 + \/3 (2) 



V3 



V3,-M (2) 

(2) 



3 i-/i-jfc. 



6 



6 



(20) 



for j ^ k. It is easy to see that \g)''s form orthogonal 
bases for subspace ^4, while |/i)'s form orthogonal bases 
for subspace S^. And there are n{n + \){n — l)/3 |.g)'s 
and n(n + l)(n — l)/3 |/i)'s. Thus we can rearrange the 
footnotes, and get \gi) and where i can change from 
1 to n(n + l)(n — l)/3. From the explicit expressions of 
\gi) and \hi), we can easily get 



Sij, 



(21) 



and and form Jordan bases for ^4 and 6*5. 

The two density operators that we want to distinguish 
can now be expressed as 



Pi 



v?-{n + 1) 



P-2 = 



7i^{n + 1) 



io 



where io = n{n + l)(n — 1). is the projection onto 
the subspace Sq, and 



^0= E i«gi)(^gii- 



(23) 



Now, let Ti be the two-dimensional space spanned by 
the nonorthogonal but linearly independent vectors \gi) 
and \hi). The Ti form a decomposition of subspace 
into n{n + l){n — 1) mutually perpendicular two- 
dimensional subspaccs. Our next problem is how to dis- 
tinguish the two pure sates in every subspace Ti. 



IV. DERIVATION AND IMPLEMENTATION OF 
THE OPTIMAL UNAMBIGUOUS 
DISCRIMINATION 

We now want to unambiguously discriminate between 
the subspaces 84^ and 5*5 in 5*6. First, we distinguish the 



Jordan Basis states \gi) and \hi) within their subspace 
Ti. In subspace Ti, the apriori probabilities of \gi) and 
\hi) are 771 and 772 respectively. Here we will use the 
method mentioned in [18]. Let us define \g^) and \h^) 
the reciprocal states of \gi) and \hi) respectively, where 



and then 



{9i\h^)^{hi\g,)^Q 



\9i) - {K\9i)\K) 

^i-mg.w 

2 1 
'.g,) + -=|/i,). 



(24) 



x/3 

\h,) 



V3' 

{9t\ht)\gt) 



^l-\{g.\hW 



9i. 



Thus, \gi) and \hi) can be rewritten as 



(25) 



\9i. 
\hi 



V3, 



(26) 



We can give a physical implementation based on Neu- 
mark's theorem [12]. Following the proposals given in 
[13-16]: (a) any pure state can be realized by a single- 
photon state and (b) following Reek's theorem [16], any 
unitary transformation can also be realized by an op- 
tical network consisting of beam-splitters, phase-shifter, 
etc. [17], we can construct an optical device, which is 
presented in Fig. la, to unambiguously discriminate be- 
tween \gi) and \hi). An additional port is initially pre- 
pared in the vaccum state |0). The operation of such a 
device is described by a unitary matrix which gives 
the probability amplitudes for a single photon entering 
via inputs [g^^) and \hi) to leave the device by outputs 
\d['^}, [D^*') and Here, the four-port optical in- 

terferometer, which is presented in Fig. lb, is used and 
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Vol' 



D 



Beam-splitter 



Detector 




Vaccum 

(a) Optical realization of u. 



□ 



Phase Shifter 



(b) Optical realization of 



FIG. 1: (a) An optical six-port interferometer used to unam- 
biguously discriminate between \gi) and \hi). (b) An optical 
four-port interferometer which is capable of realizing the two- 
dimensional unitary transformation defined by Eq. (27). 



it is capable of realizing any 2x2 unitary transformation 

U2 



sinwe*'^ coswe**^ 
coswe'^ — sinwe*^ 



(27) 



where we have tacitly assumed that the action of all the 
2x2 beam splitter upon the states of the photon is de- 
scribed by a unitary matrix of real coefficient 



sm uj cos uj 
coscj — sincj 



(28) 



The parameter u describes the transmittance {\/t = 
sinw) and the reflectivity (-y/r = cosw) of the beam split- 
ter. As shown in Fig. la, the action of this device, de- 
noted by U3 gives, 



— sinwilD^'-*) -I- coswi cos W2|£'2 
+ coswi sin 

- smuj2\D'^^) + cosu}2\F^''>) (29) 



From Eq. (26), U3 will transform l^;) and \hi) into 
V3 



■ suiuji\d[''>) 



,V3 



I /if"*) 



-I- (-^ cos wi cos 012 + - sina;2)|-D 



cos Wi smu;2 



■ cos 



— sina;2|-C'2''') 



cos LU2 



i02)\F^'^) (30) 



(31) 



As shown in Fig. la, both \gi) and \hi) have their own 
detectors Z?^''' and I?2*'' ^^'^ these two detectors will tell 
us which the input is, while F^^^ corresponds to failure. 
This suggests that no photon can be detected by detector 
^2^ (Di^) when the input is l^i) {\hi)). Thus, from Eq. 
(30) we get 



cos UJ2 



sin^ UJ2 



1 



1-1-3 cos^ uii ' 

o 2 

3 cos oji 



1 + 3 cos" cji 
Eq. (30) and Eq. (31) are reduced to 



(32) 



I ft 



V3 
2 



,V3 

-I- ( cos Ldi sm UJ2 



^^2)|F('^(33) 



3 cos ui 1^ 
1 + 3 COS uji 



1 



1 + 3 COS uji 
In other words, we can choose 



(34) 



{i)out 



n 



{i)out 



with 



nf)°"*|/i°"*> 



n 



{i)out 



n: 



{i)QUt 



{i)out, out\ 



n 



{i)out 



0, 
I, 



(35) 



(36) 



to be our detection operators. Here, also from Fig. la, 
we can get 



2 



V3 



S\IVUJ2\hj') , 



(37) 
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where Eq. (25) has been used. 

By projecting these detection operators back onto the 
space Ti, we can get 



n 



-sm^UJ2\hi){h^\ 



1 + 3 cos uji 



n'^p = I n["'' - n 



(38) 

which form the POVM detection operators for the un- 
ambiguous discrimination between \gi) and \hi). 
So we find 



(5.|nf Iff 



sin^wiKs^lg,^)!^ 



■ sm uji, 



{h^\Il'i'>\h, 



4 cos^ uji 
1 + 3 cos'^ wi 
3 cos^ wi 



1 + 3 cos wi 

And the probability of successfully identifying the two 
states is 

t(')|„ \ I ^ /u Itt(«)i 



3 2 3772 COS Wi 

-?7i sm wi H 7j — . 

4' l + 3cos2a;i 



(40) 



By letting x = 1 + 3cos^iuji(l ^ a; ^ 4), P'(wi) can be 
rewritten as 



4 a; 



(41) 



This function, (x) , has the property that dP^ (x) /dx 
happens at 



X = Xo 



(42) 



The optimal value of P*(x) denoted by p*°^* in the do- 
main 1 ^ X ^ 4, can be gotten in following three cases: 
(a) if ^ ^ 771 ^ ^ , the requirement 1 ^ x ^ 4 is satisfied, 
thus we have 



(b) if 771 < , we should choose 



P'-* . P'(x ^ 4) ^ ^,2, 



and (c) if ?7i > |, 



P 



p\x = 1) = -m. 



(43) 



(44) 



(45) 



-o 



1".)- 

|o> 
|o> 

|o>. 



U-. 



U-. 



-Oof 



Optical realization of l", in Fig. la Detector Vaccurr 



FIG. 2: Implementation of the programmable discriminator 
between pi and p2 (l^'i) and |>I'2)), which consists of each 
implementation of discrimination between \gi) and \hi). 



Our finding can be summarized as follows 



1^2 



if 771 < ^ 

1 - 2^Tnm if 5 < r/i I 



(46) 



if 7,1 > I 



Now, the discrimination between the two density op- 
erators pi and p2 is to achieve n(n + l)(7i — l)/3 dis- 
criminations described above simultaneously. But, here 
the probability for the occurrence of the input state in 

o 

Ti is piTi) = .) -. We can also give its implemen- 

n (71 + 1) 

tation in Fig. 2. We leave alone, since they are 

in both subspaces 5*1 and S2, and can't be distinguished. 
The POVM which unambiguously distinguishes Si and 



7 



S2 has the form 



Hi 



i=l 

io 



4 cos wi 
1 + 3 cos^ oji 



i=l 



Ho = I - Hi - n2, (47) 

where = n[n + l)(7i — l)/3. Thus, the probability of 
successfully identifying the two density operators pi and 
P2 is 

P(a;i) r/iTr(nipi) + ?72Tr(n2P2) 

0-2 »o 

/ Sm UJl ■r-^ I I m2 



^ ' Z — 1 



! cos^ cji 



n^(7i + 1)(1 + 3cos^ wi) ^ 
- sin^ + 2(^- l)^2COS^u;i^ ^^^^ 



2n 



+ 3 cos wi) 



Then, the optimal success probability P°p* can be ex- 
pressed as 



popt _ ^ 



n — 1 „ 



iim<l 

2(1^(1 -V^) ifl^m^l . (49) 



n — 1 



if ^71 > I 



Now, we come back to the original problem that dis- 
crimination between two pure states (l^'i)) and \'ip2) 
(|4'2)). As mentioned before, the POVM formed by IIi, 
112 and Ho with 

Hil^-a) =n2|*i) =0 (50) 

can be used. Its implementation is the same as discrimi- 
nation between pi and p2 in Fig. 2. The success proba- 
bility is 

P(wi) = 7]l(«'l|ni|«'i) +7?2(^'2|n2|*2) 

,1 .2 2772 COS^ Wi 

— in sm 0J^ + ■ 

(51) 

where the relationship 



(^7,isin^c.i+ 7:^%"^ )(i-i(7Aii7^2)n, 

1 + 3 cos uji 



:(l-|(^i|V^2)P) (52) 



1 

2* 



Input 1 1> 



|o)- 



|o).. 



-\n-l) 



|o>- 



Four-port optical 
interferometer in Fig. 1b 



|2> 



|1> 



FIG. 3: An optical device to prepare the state ]?/;) — I')- 



has been used. Finally, we can also get the optimal suc- 
cess probability as follows 



1 



popt 



' ^772(l-|(V'l|V^2)P) if^l<5 

|(l-Vm^)(l-|(V^l|^2)P) if ^^771*^1 , 
. ^7;i(l- 1(^1 17^2) P) 



which is apparently independent of n, the dimension of 
space H, as we mentioned at the beginning. 



V. DISCUSSION AND CONCLUSION 

How to prepare the states with single photon is an- 
other important question in the optical realization of the 
discriminator. We shall give a brief discussion about this. 
An optical setting shown in Fig. 3 can prepare the state 
1^) = S"=i OiK)- It is constructed by a series of 2x 2 uni- 
tary transformations which can be realized by the four- 
port optical interferometer in Fig. lb. If the parameters 
of each four-port interferometer are fixed properly, this 
setting achieves a single-photon state in n-dimensional 
space H. We denote the operation of this device by 
and 



U\i>)\i) 



(54) 
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Then, by using U\^^^, U\^,,^ and (7|^,) in succession, we 
can get 

|*?> = IV'i) ® IV'?) ® IV'2) 

= C/|^2)?7|^,)C/|^^)|1) 

= C/|^^)(»?7|^,)®?7|,0,)-|1)®|1)®|1) (55) 

where lip'?) is either \ipi) or IV'2) and |1) (g) |1) (g) |1) denote 
the only input for the device. And j^"!) or \'i>2) is to be 
prepared depending on ?7|^,) being C/|^i) or U\^,.^y 

In conclusion, we have reconsidered the problem of 
the universal programmable quantum state discrimina- 
tor originally introduced in [1]. And we solved a more 
generalized problem that the unknown states are qudit 
states in the n-dimensional {n ^ 2) Hilbert space H. 



We adopted the equivalence between the discrimina- 
tion of unknown pure quantum states and that of known 
mixed states. With the Jordan-basis method, we simpli- 
fied the problem and finally achieved the optimal unam- 
biguous discrimination between the two unknown states, 
and the corresponding detection operators have already 
been given. Significantly, we arrived at an important 
conclusion that the optimal success probability of the 
discrimination between two unknown pure states is inde- 
pendent of n, the dimension of the space H. 

Furthermore, we also give the implementation of the 
optimal POVM based on the Neumark's theorem. The 
POVM can be implemented on a larger Hilbert space, 
where the additional degrees of freedom called ancilla 
are needed. 
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